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ABSTRACT 

We describe detailed simulations of X-ray-emitting populations to evaluate 
the levels of contamination by both Galactic and extragalactic X-ray sources un- 
related to a star-forming region under study. For Galactic contaminations, we 
consider contribution from main-sequence stars and giants (not including cat- 
aclysmic variables and other classes of accretion-driven X-ray binary systems) 
as they make the dominant contribution at the position of the Carina Nebula. 
The simulations take into consideration a variety of technical factors involving a 
Galactic population synthesis model, stellar X-ray luminosity functions, Chandra 
telescope response, source detection methodology, and possible spatial variations 
in the X-ray background and absorption through molecular clouds. When ap- 
plied to the 1.42 square-degree field of the Chandra Carina Complex Project 
(CCCP), the simulations predict ~ 5000 contaminating sources (1 source per 
square arcminute of the survey), evenly distributed across the field. The results 
of the simulations are further employed in a companion CCCP study to assign 
membership probabilities to individual sources. 
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Subject headings: Galaxy: disk; ISM: individual (Carina) - open clusters and 
associations: individual (Carina) - stars: formation - stars: pre-main sequence - 
X-Rays: stars 



1. Introduction 



Identifying stellar members associated with a star-forming region is important for studies 
of the stellar initial mass function (IMF), cluster dynamics, and star formation processes. 
Stars with dusty protoplanetary disks are readily found in infrared (IR) surveys where disk 
emission is strong, but when the disks have dissipated, the stars are most easily identified 
by their strong X-ray flaring. High spatial resolution X-ray observations of an active star 
forming region with the Chandra X-ray Observatory wil l often reveal hundred s or thousands 
of pre-main sequence (PMS) stars, as well as OB stars ( iFeigelson et al.l 120071 ). 



An important challenge for membership studies is to distinguish young stars formed in 
a molecular cloud from various types of contaminant sources along the line-of-sight. In X- 
ray survey s of star-forming regions, stellar X -ray contaminants include main -sequence (MS) 
stars (e.g. ISchmitt et al.lll995t ISchmittI 119971 ) and some types of giants (e.g. iPizzolato et al. 



2000; Gondoin 



2005! ) . These foreground and background stellar contaminants have much 



less impact on surveys in X-ray images than IR images b ecause magnetic activity in PMS 



stars is elevated 10^ — 10^ above levels in older stars (e.g. iPreibisch &: Feigelsonll2005l ). For 
star forming regions located in the quadrant of the Galactic plane centered on the Galactic 
Cent er, contamination by cataclysmic variables should be also considered (e.g. iMuno et al. 



2004 ). X-ray surveys also suffer some contami nation by extragalactic so urces, mainly quasars 
and other active galactic nuclei (AGN) (e.g. iBrandt fc Hasingerll2005l ). These can be seen 
through the Galactic Plane as faint, absorbed X-ray sources. 

Careful simulations of contaminating populations can provide estimates for the number 
of contaminants, their spatial distribution, and some of their apparent X-ray and IR prop- 
erties. T hese estimates can inform efforts to assign individual classifications to the detected 
sources ( iBroos et al.ll2011b[ ). Such simulation studies have been performed for the Orion 



Nebu l a Clu ster (ONC) and the Cepheus B (Cep B) star-forming regions (iGetman et al. 
20051 . booeh . For the deep Chandra observation of the ONC, 150 - 200 AGNs, 15 - 20 
foreground stars, and zero background stars are predicted, constituting ~ 10 — 15% of the 
> 1600 X-ray detected sources. For the shallow Chandra observation of Cep B, 20 — 30 
AGNs, 10 — 15 foreground stars, and a few background stars are predicted, again ~ 10% of 
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the Chandra source population. 



The number and observed properties of contaminating populations depend on a num- 
ber of factors: observatory sensitivity, source detection method, observation exposure time, 
survey field of view. Galactic direction, distance to the star forming region, and absorption 
through the local molecular cloud. Therefore, it is ill-advised simply to extrapolate the con- 
tamination results for the Cep B and ONC Chandra fields to substantially different X-ray 
observations. 

The purpose of this study is to present our contaminant simulation method a nd its 

The 



application to the Chandra Carina Complex Project fC CCP: iTownslev et a. 



simulation methods are refined from those described by iGetman et al.l (|2006[ ). The Carina 



2011ai ) 



Nebula (NGC 3372) is one of the richest nearby high-mass star-forming co mplexes, located 
2.3 kp c from the Sun towards the Galactic direction (/, b) = (287°. 6, — 0°.8) (ISmith fc Brooks 
20081 ). The CCCP combines 22 spatially contiguous C handra-AClS-l pointings with nominal 
exposure of 60 ksec coverin g 1.42 square degrees a round the chain of famous massive star- 
forming regions in Ca rina (ITownsley et al.l |2011aJ) . A relatively shallow molecular cloud 
(lYonekura et al.lboosi ) covers roughly 1/3 of the CCCP field (Figure [It). More than 14000 
CCCP point sources have been detected (IBroos et al.ll2011al ). 



The simulations of Galactic stellar and extragalactic contaminating populations in 
CCCP are described in §5J They show that ~ 5000 contaminating sources are expected 
in the CCCP catalog. Basic properties of these simulated contaminating sources (spatial 
distribution. X-ray median energies, and J-band magnitudes) are presented in ^ The 
comparison of the properties of the simulated contaminating sources and the real CCCP 
sources is considered, at a simplistic level, in §11 Based on the results of this study, with 
additional information on X-ray variability, IR excess, and optic al spectral types, ind ividual 
membership classifications for the CCCP sources are derived by iBroos et al.l ( 1201 Ibl ). 



2. Simulations of Cluster Contaminants 



2.1. Extragalactic Sources 



We evaluate the ex pected contamination by ex tragalactic X-ray sources using methods 
similar to the analyses of iGetman et al.l ( I2005l . l2006l ). Nominal (corresponding to the the high 
Galactic latitude fields with low Galactic absorbing column used for extragalactic studies) 
hard-band (2 — 8 keV) fluxes incident on Chandra are drawn fr om the log — log S distribu- 
tion of extragalactic sources described by iMoretti et al.l (l2003l) assuming a pow er law source 
spectrum consistent with fiux-dependencies described by iBrandt et al.l (l200l[ ). More than 
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22600 extragalactic sources down to a nominal flux of \og{Fx) = —16.7 erg cm~^ s~^ are pre- 
dicted for the 1.42 deg^ CCCP field. One hundred Monte Carlo simulations are constructed 
by placing > 22600 artificial sources randomly across the CCCP field. Source photon indices 
are drawn from uniform distributions in the ranges —0.5 < F < 2 and 1 < F < 2 for nominal 
fiuxes \og{Fx) < —14.1 erg cm~^ s~^ and \og{Fx) > —14.1 erg cm~^ s~^, respectively. 

The nominal fiux of each simulated source is transformed to apparent (i.e., observed) 
fiux, taking into consideration an extra absorbing column that is estimated as the sum of 
two components: a uniform HI column de nsity of Nfj ~ 1.4 x 10^^ cm~^ through the entire 



Galactic disk in the direction of Carina ( Dickey fc Lockman 1990 1^ and spatially variable 
absorption through the Carina molecular clouds. The cloud component is derived from 
th e velocity-integrated ra aps of CO line emission obtained with the NANTEN telescope 
by lYonekura et al.l (120051 ). The gas column density outside the dense molecular cores was 
estimated from the maps of ^^CO emission intensity integrated over the velocity range of 
30 to - lOkms"^ (K ^^CO) [K km s'^]) using the relationsh ip N(Hy) = I ( ^^CO ) x 1.6 x 10^° 



(lYonekura et al.l l2005l ) . Adopting the gas-to-dust ratio of IVuong et al.l (120031 ) , these gas 
column densities correspond to visual extinctions up to Ay ~ 3 mag in the molecular clouds. 
The dense core column densities are evaluated from th e C^^O measurements and correspond 
to Ay ~ 5 — 10 mag (Table 3 in lYonekura et al.ll2005l ). 



Chandra source photon spectra were then simulated, and observed count rates and 
median energi es were derived using the fakeit command of the XSPEC software package 
(lArnaudlll996l ). This program samples photons from the AGN spectrum, calculating the soft 
X-ray absorption by the line-of-sight material, and convolving the results with the Chandra 
mirror and ACIS detector response. Typical on-axi s CCCP calibration redistribution matrix 
and auxiliary response calibration files were used (IBroos et al.l 1201 lal ). After applying local 
background levels found in the CCCP data, we then remove very weak extragalactic sources 
that would have fallen below the CCCP source detection thresholcj^. 



^This was obtained from NASA's HEASARC tool located at http: //heasarc .gsfc .nasa.gov/cgi-bin/Tools/w3nh/w3n! 

^Position-dependent detection thresholds corresponding to the CCCP detection process were calcu- 
lated using local background levels estimated from t he CCCP data. The ver y sophisticated and in- 



novative source detection technique used in CCCP (jBroos et all l201Cll . l2011al) can be approximated 
by thresholding the source significance statistic, Signif, which is calculated by the ACIS Extract 
(http://www.astro.psu.edU/xray/docs/TARA/ae_users_guid.e.html) package as the ratio of net counts 
to the uncertainty on that quantity. The distribution of source significance for all CCCP sources peaks at 
Signif = 0.8, has a very sharp decay below 0.8 encompassing only 3% of sources, and has an exponen- 
tial decay above 0.8. In our contamination simulations we thus choose the source detection threshold of 
Signif — 0.8. 
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Our simulations predict that ~ 2500 extragalactic sources should be detected in the 
CCCP observation. 



2.2. Galactic Sources 



It has been argued that the point source component of the hard "Galactic ridge emis- 
sion", which stretches over 90 degrees in Galactic longitude and only 2 degr ees in latitude 
around the Galactic Center, is mainly compose d of cataclysmic variables (e.g lEbisawa et al. 
20051 : iRevnivtsev et al.ll2007l : iHong et al.ll2009[ ). The current version of our simulation omits 
the distribution of hard (> 3 keV) X-ray sources in the Galactic plane attributed to cata- 
clysmic variables and other classes of accretion-driven X-ray binary systems. However, this 
population and its X-ray emission should be unimportant in the Carina star forming region 
because it is located 72 degrees away from the Galacti c Center where the Galactic Ridge 
emission is weak (see Figure 1 in IRevnivtsev et al.ll2006l ). 



We simulate the stellar contamination from main sequence and giant stars in the Galac- 
tic Plane in three stages. First, Monte Carlo simulations of the Galactic stellar population 
expected wit hin the CCCP fie ld are constructed based on the 'Besangon model' of stellar 
populations ( iRobin et al.ll2003l ). The simulations are computed by the Web service provided 
by the Besangon group|f|. Using models of Galactic structure, stellar mass function, and 
star formation history, the calculation gives the number distribution of stars along a cho- 
sen line-of-sight including individual spectral types, optical and near-infrared magnitudes, 
ages, distances, and source extinctions. An additional position-dependent absorption rep- 
resenting the Carina molecular clouds was added for stars behind Carina {D > 2.3 kpc). 
Besangon model limitations include a simplistic assumption of uniform extinction through- 
out the Galactic Plane and omission of spiral arms. We run the Galactic population model 
with an unrealistically deep V = 28 mag detection limit, as the faint stars are later omitted 
by X-ray selection. 

Second, X-ray luminosities for each star obtained from a Besangon model simula- 
tion are then estimated by Monte Carlo sampling of X-ray luminosi ty functions (XLFs) 



of MS stars measured in the solar neighborhood from i^OS'^T surveys (jSchmitt et al.lll995 



SchmittI 119971 : iHiinsch et al.l Il999l ). Based on these studies, we adopt the following log- 



normal XLFs with the following means and standard deviations in logLx [erg s ^] in 
the ROSAT PSPC energy band (0.1 - 2.4 keV): < 28.20 > ±0.79 for MS F-type stars; 
< 27.94 > ±1.02 for MS G-type stars; < 27.86 > ±0.65 for MS K-type stars; and 



• ^http : //model . obs-besancon . f r/[ 
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< 27.44 > ±0 76 fo r MS M-type stars. For giant stars, we adopt < 28.91 > ±0.86 
( Pizzolato et al. X-ray spect r al sha pes were assumed to follow the stellar temperature- 



luminosity relations of iGiidel et al.l (Il998[ ) and iPizzolato et al.l (120001 ) for MS and giant stars 



respectively. The adopted temperature-luminosity relations in the ROSAT (0.1 — 2.4 keV) 
band are: log(T) [MK] = -9.2 + 0.331og(Lx) [erg s"^] for MS stars and log(r) [MK] = 
—5.1 + 0.21og(Lx) [erg s""*^] for giants. 

Third, given the X-ray luminosity in the ROSAT (0.1 — 2.4 keV) band, spectrum, dis- 
tance, and line-of-sight absorption for each of the simulated stars, the count rate in the 
Chandra ACIS detector in the (0.5 — 8.0 keV) band is obtained using the Portable In- 
teractive Multi-Mission Simulator (PIMM^Ifl The detectability of each simulated star is 
then evaluated by locally-defined thresholds of the source significance statistic, similar to 
our treatment of extragalactic contaminants (§ 12.11) . For the stars that satisfy the detection 
threshold, the XSPEC fakeit command was used to simulate the spectrum to derive X-ray 
median energies for individual stars. 

A typical run of the Besangon Galactic disk population model of the 1.42 deg^ CCCP 
exposure of the Carina Nebula region predicts ~ 200, 000 foreground (distance < 2.3 kpc) 
MS stars in the following mass distribution: 4,000 F stars, 10,000 G stars, 20,000 K stars, and 
170,000 M stars. About 1,000 foreground giants are present. Behind the cloud at distances 
2.3 - 5.0 kpc, the model predicts ~ 930, 000 background stars: 30,000 F stars, 77,000 G 
stars, 180,000 K stars, 630,000 M stars, and 9,000 giants. The vast majority of these stars 
are undetectable in the CCCP X-ray exposures. This is expected from the appearance of star 
forming regions at distances around 1 — 3 kpc in different wavebands: infrared images are 
typically dominated by dense populations of field stars, while X-ray images are dominated 
by PMS stars in the young clusters. 

A typical Monte Carlo run predicts ~ 1800 foreground and ~ 900 background field stars 
will be detected in the X-ray survey. The foreground detections include: 150 F stars, 450 
G stars, 150 K stars, 850 M stars, and 150 giants. Approximately 20%, 50%, and 30% of 
the simulated foreground detections have distances from the Sun in the ranges < 0.5 kpc, 
0.5 — 1.5 kpc, and 1.5 — 2.3 kpc, respectively. The background detections include: 100 F stars, 
500 G stars, 10 K stars, 40 M stars, and 250 giants. Approximately 35%, 35%, and 30% of 
the simulated background detections have distances in the ranges 2.3 — 3.0 kpc, 3.0 — 4.0 kpc, 
and 4.0 — 5.0 kpc, respectively. 



^http : //cxc . harvard . edu/ciao/eihelp/pimins . html 

^At this stage of the simulations, it was not feasible to use XSPEC to derive count rates for our very 
large sample of simulated stars. 
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Several caveats pertain to the simulation. First, the hne of sight towards the Carina com- 
plex is tangent to the Sagittarius-Carina spiral arm, but the Besangon model of the Galactic 



disk population (which our simulations rely on) lacks spiral arm modeling (IRobin et al. 
20031 ). Second, X-ray luminosities of stars with ages < 1 Gyr might be even higher than 



those of stars from the solar neighborhood (jPreibisch fc FeigelsonI 120051 ). Roughly 1/4 of 
the simulated Besangon stars have ages < 1 Gyr. Third, due to computational reasons, the 
simulations of Galactic background stars were truncated at a distance of 5 kpc. Due to 
these three factors, we might expect even larger stellar contamination populations. Fourth, 
on the other hand, a possible decrease in X-ray activity with mass and age in giant stars 
( iPizzolato et al.ll2000[ ) is ignored in our simulations. If one considers only the most plausible 
class of X-ray emitting giants, F and G giants with ages < 1 Gyr, the expected number of 
detected giants in the CCCP field could be as small as ~10 in the foreground and ~10 in 
the background. 

Finally, we note that many Galactic field stars will each produce < 3 photons in the 
CCCP image, too faint for detection as a point source. The undetected foreground stars 
produce a total of ~ 10, 000 X-ray photons in the soft (0.5 — 2 keV) band. These photons 
will appear as a uniform diffuse X-ray component in the CCCP image. These stars likely 
contr ibute ^1% of the obse rved > 10^ count soft X-ray emission that pervades the CCCP 
field (ITownsley et al.ll2011b[ ). 



3. Properties of Cluster Contaminants 



Basic properties of the simulated contaminating populations — spatial distributions. 
X-ray mediar i energies, and J-ba. nd magnitudes — are described here. These properties are 
employed by iBroos et al.l (l2011b( ) to establish Carina membership probabilities for each of 
the CCCP sources. 



3.1. Spatial Density and Distribution 



Our simulations predict that roughly 5000 contaminants, or approximately 1 source 
arcmin"^ in the 5112 arcmin^ CCCP survey, will be detected in our CCCP observations. 
About half of these are expecte d to be extragalact ic sources, and half Galactic field stars. 
With > 14000 detected sources (IBroos et al.ll2011al ). contaminants are expected to comprise 
~ 35% of the detected CCCP source population. For comparison, the source surface densities 
of the X-ray contaminants detected in the ONC and Cep B observations (§ [T]) are 0.8 and 



- 8 - 



0.2 source arcmin"^, respectively. A variation among these numbers is a product of multiple 
competing factors such as, detection technique, exposure times, thicknesses of molecular 
clouds and distances to the clouds. Galactic coordinates. 

Across most of the Carina molecular cloud, the absorpt ion is low, typically A y < 3 mag, 
and reaches ~10 mag only in a few small molecular cores ( lYonekura et al.l 120051 ). Thus the 
spatial distributions of the simulated contaminants for all classes of the contaminants are 
relatively uniform across the field of view. We ignore small differences in the surface density 
of field stars between the northern (closer to the Galactic plane) and the southern portions 
of the CCCP field, predicted by the Besangon model: < 2% and < 5% for the foreground 
and background stars, respectively. The spatial distribution of the contaminants is very 
different from the clustere d distribution of the observed CCCP source s, as shown in maps 
of the source distributions fiTownsley et al.ll2011ai iFeigelson et al.ll201ll ). The difference can 
be quantitatively seen in projection along the right ascension axis (Figure [1]). The lower 
numbers of simulated contaminating sources at the both ends of the right ascension range, 
compared to that of the central part of the field (Figure [1^), are simply due to the decreasing 
angular coverage of the CCCP observation near the east and west boundaries (Figured}:;). 



3.2. X-ray Median Energy 

The X-ray median energy statistic, MedE, is the most robust among the measured 
photometric quantities used to describe the spectral shapes of weak X-ray sources. MedE 
is most effective as a surrogate for the absorption column, but if the absorption is indepen- 



dently known, it ca n give a rough measure of plasma temperature (IFeigelson et al.l 12005 



Getman et al.ll2010l ). From the simulated spectra of the Carina contaminants (§ [2]), we ob- 
tain histograms of MedE for each of the three simulated contaminant classes and compare 
them with the observed CCCP sources in Figure [2l Notice that the observed CCCP sources 
are mixture of objects from all four classes: PMS stars, foreground stars, background stars, 
and AGNs. A few important findings emerge. 



1. The three contaminant classes occupy mostly distinct ranges of MedE: 96% of the 
simulated foreground stars have very soft X-ray spectra with 0.7 < MedE < 1.1 keV 
range; 81% of the simulated background stars are within the 1.1 < MedE < 1.5 keV 
range; 94% of the simulated extragalactic sources are within the 2.5 < MedE < 4.5 keV 
range. 

2. The observed CCCP sources span a wide range of MedE, from 0.6 to 7 keV, with a 
skewed distribution peaking at 1.4 — 1.5 keV. 
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3. The histograms of the observed CCCP sources and the simulated foreground stars have 
similar numbers of sources in the range 0.7 < MedE < 1 keV, suggesting that nearly 
all the CCCP sources in this range are likely foreground stars. This also implies that 
our simulation did not significantly underestimate the foreground population due to 
an unmodeled population of MS stars in the Carina spiral arm (§ 12. 2p . 

4. The histograms of the observed CCCP sources and the simulated AGNs have similar 
numbers of sources in the range 2.5 < MedE < 4.5 keV, suggesting that nearly all the 
CCCP sources in this energy range can be AGNs. 

5. The peak of the simulated background stars coincides with the peak of the observed 
CCCP sources, but the numbers of the CCCP and simulated background sources are 
drastically different. It is clearly impossible to distinguish between background stars 
and Carina PMS members using MedE. 

We see that the X-ray median energy is an important discriminant for distinguishing 
some contaminant populations. It can be effective for foreground stars and extragalactic 
sources, but not background stars which have the same median energies as many Carina 
PMS stars. 



3.3. J-band Magnitude 



Younger PMS stars with protoplanetary disks (infrared Class O-I-II systems) can be 
readily distinguished from Galactic field stars by their infrared excess. But often, as in the 
Carina complex, the majority of stars associated with a star forming region are older, diskless 
(Class III) systems. X-ray studies are particularly effective in detecting Class III stars which 
have high levels of magnetic activity. While in some X-ray studies of star forming regions, 
IR color-color and/or color- magnitude dia grams can be employed to successfully distin guish 
between diskless PMS and older stars (e.g. iGetman et al.ll2005l . l2006l : iKuhn et al.ll2010l ). the 
small source extinction range (typically only Ay ~ 1 — 3 mag for most parts of the Carina 
region) and the large range of distances (only with respect to color-magnitude diagrams) for 
the simulated contaminating stars ( §2.2p blunt this technique. Therefore, here we restrict our 
consideration to a simple near - IR pho tometric property, J-band magnitude. In addition to J- 
band magnitude, iBroos et al.l ( l2011bl ) employ some near-IR and mid-IR colors to distinguish 
between disk-bearing PMS and other stars. Wider usage of IR colors and incorporation of 
IR-to-X-ray flux ratios (when applicable) will be considered in our future studies of other 
star forming regions. 
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Figure [3] compares J-band magnitude histograms among the 6, 182 CCCP sources with 
available 2MASS counterparts within the entire CCCP field of view (dashed lines) and the 
three simulated contamination populations (solid lines). J-band magnitudes for the fore- 
ground and background field stars are obtained directly from the Besangon model simula- 
tions; for the background stars, magnitudes are slightly increased for the absorption from 
the Carina molecular cloud. For the extragalactic contaminants, we produce an approxi- 
mate distribution of J-band magnitudes using the relationship between the nominal (see 
§2.ip hard-band X-ray fiu x and /-band rnagnitu de reported for the Chandra Deep Field 
North sources (Figure 7 in [Alexander et al.ll200ll ). The /-band magnitudes are transformed 
to J-band magnitudes following Figure 11 in Alexander et al., accounting for an average 
combined visual absorption of Ay ~ 10 — 12 mag through the Galactic disk in the direc- 
tion of Carina and through the Carina molecular clouds. For better compatibility with 
the CCCP histograms, we also present (Figure (Sb) the histograms of J-band magnitudes 
for simulated foreground and back ground stars mu l tiplie d by the "generic" 2MASS J-band 
detection completeness curve from Skrutskie et al. ( 2006 1^. 



Several results emerge from the J-band magnitude distributions. First, the histogram 
of the simulated foreground population has two peaks (red curve). The fainter peak from 
M-type field dwarfs, constituting ~ 40% of the simulated foreground population, lies below 
the 2MASS sensitivity limit of J ~ 16.5 — 17 mag. A bimodal shape is indeed seen in 
the deep er HAWK-I J-band m agnitude histogram of the CCCP sources with MedE < 
1.1 keV (jPreibisch et al.l l201ll ). Second, for bright stars with J ^ 16 mag, the J-band 
magnitude histogram of the simulated foreground population (red curve) agrees well both in 
number and shape with the magnitude distribution of CCCP sources with median energies < 
1.1 keV (black dashed curve). This indicates, in agreement with the \ow-MedE distributions 
discussed in § 13.21 that most of the observed CCCP sources in this sub-sample can be 
attributed to foreground stars. Third, roughly 60% of the simulated stellar background 
population lies below the 2MASS sensitivity limit of J > 16.5 mag. Fourth, most of the 
simulated extragalactic objects lie below J = 21 mag. The HAWK-I survey can detect some 
of these extragalactic objects near its detection limit (jPreibisch et al.ll201ll ). 



^See also the 2MASS J-band detection completeness curves in the 2MASS All Sky Data Release Documen- 
tation at http : //www . ipac . caltech . edu/2mass/releases/allsky/doc/sec6_5al . html Notice, however, 
that this "generic" 2MASS detection completeness curve does not account for the effects specific to the CCCP 
field, such as the source crowding and IR nebulosity. 
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4. Spatial Distributions Stratified by Median Energy and J-band Magnitude 



Here we consider the spatial distribution of CCCP sources stratified by MedE and by 
J-band magnitude that, as seen in § [3l are often linked to contaminant classes. We expect 
that Galactic and extragalactic contaminants will have random locations in the CCCP field, 
but many Carina members will lie in young stellar clusters. The spatial distributions are 
shown in Figures HHSj with MedE strata informed by Figure [21 The contours sh ow the 
distribution of molecular material from NANTEN CO maps (lYonekura et al.l l2005l ) . The 
CCCP sources are further divided into two groups that have or lack 2MASS counterparts, 
giving two mag nitude strata divided approximately at J ~ 16.5 mafl Interpretattan ot 
spatial distributions of CCCP sources is co nfused by the large variations in sensitivity due 
to overlapping exposures and off- axis angle (IBroos et al.ll2011al ). Our qualitative discussion 
here is designed to assess the role of contaminants. A more quantitative trea tment of spatial 



group ings from a spatially complete sub-sample of CCCP sources is given by lFeigelson et al. 
fl201lh . 



4.1. Hard CCCP Sources 

In this section we consid er all CCCP sources with X-ray median energy MedE > 2 keV, 



equivalent to Ay > 10 mag (iGetman et al.l 120101 ). In panels d — f oi Figure HI 70% of the 



~ 700 CCCP sources with 2MASS counterparts are spatially concentrated towards the edges 
of the Carina molecular clouds. Three prominent source groupings are seen: members of the 
Treasure Chest cluster, stars within the Trumpler 16 cluster, and stars at the interface 
between the rich Trumpler 14 cluster and C^^O core ^^10. The latter region also has a dense 
grouping of MSX mid-infrared sources. These groupings of absorbed X-ray sources are likely 
to be very young stellar objects. In contrast, out of ~ 2800 hard CCCP sources without 
2MASS counterparts (Figure — /), ~ 2000 are evenly distributed across the CCCP field 
and are likely extragalactic candidates. This number is in satisfactory agreement with the 
number of simulated AGNs (2500 from §2.11) . 



^ A ~ 3' circular 'hole' in the distribution of 2MASS sources is seen in Figure due to saturation by rj 
Car in the 2MASS images. 
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4.2. Soft CCCP Sources 

In this section we consider all CCCP sources with X-ray median energy MedE < 2 keV. 
Figure m presents all CCCP sources with available 2MASS counterparts. Panel a {MedE < 
1.1 keV) shows that only 40% of the sources are spatially concentrated in the ridge of fa- 
mous young stellar clusters in C arina (including Trum p ler 15, Trumpler 14, Trumpler 16, 



Treasure Chest, and Bochum 11 (IFeigelson et al.ll201ll ) rl. while the majority of the sources 



are relatively uniformly distributed across the CCCP field of vie'w^ (ignoring the "egg-crate 



effect", due to variations in sensitivity with off- axis angle (IBroos et al.l l2011al )). In con- 
trast to panel a, 70% of the sources in panel b (1.1 < MedE < 1.5 keV) and 80% of the 
sources in panel c (1.5 < MedE < 2.0 keV) are seen in projection against the chain of 
Carina clusters. Figured presents all CCCP sources without 2MASS counterparts. Panels a 
{MedE < 1.1 keV), b (1.1 < MedE < 1.5 keV), and c (1.5 < MedE < 2.0 keV) show that 
approximately 80%, 50%, and 50% of sources, respectively, are evenly distributed across the 
CCCP field. 

Thus, based simply on their uniform spatial distributions, as many as 700 (+1400) 
with MedE < 1.1 keV, 800 (+1100) with 1.1 < MedE < 1.5 keV, and 300 (+700) with 
1.5 < MedE < 2.0 keV sources with (without) 2MASS counterparts could be potential 
Galactic stellar contaminants. Further consideration of the expected MedE and J-band 
properties for the simulated stellar contaminants (§2]) suggests that indeed, most of the 
700 + 1400 = 2100 sources with MedE < 1.1 keV could be foreground stars, and roughly 
half of the 800+1100 = 1900 sources with 1.1 < MedE < 1.5 keV could be background stars. 
Therefore, this leaves us with at least ~ 300 {MedE < 1.1 keV), ~ 1000 (1.1 < MedE < 
1.5 keV), and ~ 1000 (1.5 < MedE < 2.0 keV) CCCP sources evenly distributed across 
the field that could be either additional contaminants unaccounted for by the simulations, 
and/or young stellar objects found outside the famous Carina clusters. In the former case, 
our underestimation of the Galactic stellar contamination population could be due to un- 
modeled young MS stars associated with Galactic spiral arm toward Carina, or background 
field stars with distances > 5 kpc ( §2.2p . In the latter case, we have an indication for the 
presence of a widely distributed population of young stars in the Carina complex. 



^The fie ld of the cluster ridge is defined here as the lowest contour of the X-ray source density shown in 
Figure 1 of IFeigelson et al.l (|2011l ). 



^In this section, source density estimates of evenly distributed population s are based on average source 
densities in the four ACIS-I fields outside the chain of the famous clusters. In lTownslev et al.l (|2011al ) these 
four fields are labelled as E2, E4, SBl, and SB2. 
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5. Summary 



This work describes simulations of X-ray-emitting source populations that will appear 
in high-resolution Chandra X-ray Observatory studies along the Galactic Plane. The effort 
is particularly relevant to discriminating PMS stars in star forming regions from Galactic 
field star and extragalactic contamination. Our simulations account for two X-ray classes: 
extragalactic sources (primarily AGNs) seen through the Galactic disk and Galactic field 
stars (main sequence and giants from types F to M) distributed throughout the disk. The 
simulations are based on a variety of technical considerations (§ |2]): use of a Galactic pop- 
ulation synthesis model and measured X-ray luminosity functions to obtain field star X-ray 
source distributions; convolution of realistic contaminant X-ray intensity and spectral distri- 
butions through the telescope response; consideration of spatially varying absorption across 
the field due to molecular clouds; application of a realistic source detection method including 
threshold effects due to position-dependent point spread functions and X-ray background. 
The current version of the simulations omits contributions from accretion binary systems 
and will thereby underestimate the population at high MedE for fields close to the Galactic 
Center. 

The results o f the simulations are ap plied to the CCCP X-ray survey of the Carina star 
forming complex ( iTownsley et al.ll2011al ) to evaluate the levels of non-Carina contaminants 
among the > 14000 CCCP point sources. The simulations predict that about 2500 AGNs 
( §2.ip . 1800 foreground stars, and 900 background stars ( §2.2p will be detected. The simu- 
lated contaminants are expected to have uniform spatial distributions and exhibit significant 
differences in their distributions of the X-ray median energy and J-band magnitude (§ |3]). 
These properties are thus useful for classifyin g individual CCCP sources as likely Carina 
members or contaminants ( iBroos et al.ll2011bl ). 



Comparing properties of the simulated contaminating sources and the CCCP sources, 
at a simplistic level we find that the number of hard CCCP sources evenly distributed 
across the field is in satisfactory agreement with our extragalactic contamination prediction, 
but the number of evenly-distributed soft CCCP sources is a factor of 2 higher than our 
Galactic contamination prediction (SI)- This discrepancy may be due either to unmodeled 
effects (e.g., stars beyond D = 5 kpc, an additional population of MS stars in the Carina 
spiral arm, or evolution of X-ray activity in older stars), or the presence of a widely dis- 
tributed population of young stars in the Carina complex. The inference of distributed star 
formation, or the kinematic drifting of an older generation of Carina stars, i s substantiated 



by rnore detailed studies of the CCCP sources ly ing outside the rich clusters ( iFeigelson et al. 
201 ll : IPreibisch et"aDl201ll : IPovich et allboilbh . 



- 14 - 



We thank the anonymous referee for his time and many useful comments that improved 
this work. This work is supported by Chandra GO grant SAO G08-9131X (L. Townsley, 
PI) and the Chandra ACIS Team contract SV4- 74018 (G. Garmire, PI), issued by the 
Chandra X-ray Center operated by the Smithsonian Astrophysical Observatory for and on 
behalf of NASA under contract NAS8-03060. M.S. P. is supported by an NSF Astronomy and 
Astrophysics Postdoctoral Fellowship under award AST-0901646. This publication makes 
use of data products from the Two Micron All Sky Survey (a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Center/California Institute of 
Technology, funded by NASA and NSF). 

CXO (ACIS) 

REFERENCES 

Alexander, D. M., Brandt, W. N., Hornschemeier, A. E., Garmire, G. P., Schneider, D. P., 
Bauer, F. E., & Griffiths, R. E. 2001, AJ, 122, 2156 

Arnaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis Software and 
Systems V, ed. G. Jacoby & J. Barnes (San Francisco, CA: ASP), 17 

Brandt, W. N., et al. 2001, AJ, 122, 2810 

Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 827 

Broos, P. S., Townsley, L. K., Feigelson, E. D., Getman, K. V., Bauer, F. E., & Garmire, 
G. P. 2010a, ApJ, 714, 1582 

Broos, P. S., et al. 2011, ApJS, submitted (CCCP Catalog Paper) 

Broos, P. S., et al. 2011, ApJS, submitted (CCCP Classifier Paper) 

Dickey J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Ebisawa, K., et al. 2005, ApJ, 635, 214 

Feigelson, E. D., et al. 2005, ApJS, 160, 379 

Feigelson, E., Townsley, L., Giidel, M., & Stassun, K. 2007, in Protostars and Planets V, ed. 
B. Reipurth, D. Jewitt, & K. Keil (Tucson, AZ: Univ. Arizona Press), 313 

Feigelson, E. D., et al. 2011, ApJS, submitted (CCCP Clustering Paper) 



- 15 - 



Getman, K. V., Feigelson, E. D., Grosso, N., McCaughrean, M. J., Micela, G., Broos, P., 
Garmire, G., & Townsley, L. 2005, ApJS, 160, 353 

Getman, K. V., Feigelson, E. D., Townsley, L., Broos, P., Garmire, G., & Tsujimoto, M. 
2006, ApJS, 163, 306 

Getman, K. V., Feigelson, E. D., Broos, P. S., Townsley L. K., & Garmire, G. P. 2010, ApJ, 
708, 1760 

Gondoin, P. 2005, A&A, 444, 531 

Giidel, M., Guinan, E. F., & Skinner, S. L. 1998, in ASP Conf. Ser. 154, Cool Stars, Stellar 
Systems and the Sun, ed. R. A. Donahue & J. A. Bookbinder (San Francisco, CA: 
ASP), 1041 

Hong, J. S., van den Berg, M., Grindlay J. E., & Laycock, S. 2009, ApJ, 706, 223 

Hiinsch, M., Schmitt, J. H. M. M., Sterzik, M. F., & Voges, W. 1999, A&AS, 135, 319 

Kuhn, M. A., Getman, K. V., Feigelson, E. D., Reipurth, B., Rodney, S. A., & Garmire, 
G. P. 2010, ApJ, 725, 2485 

Moretti, A., Campana, S., Lazzati, D., & Taghaferri, G. 2003, ApJ, 588, 696 

Muno, M. P., et al. 2004, ApJ, 613, 1179 

Pizzolato, N., Maggio, A., & Sciortino, S. 2000, A&A, 361, 614 
Povich, M. S., et al. 2011, ApJS, submitted (CCCP IR YSOs Paper) 
Preibisch, T., & Feigelson, E. D. 2005, ApJS, 160, 390 
Preibisch, T., et al. 2011, ApJS, submitted (CCCP HAWK-I Paper) 

Revnivtsev, M., Sazonov, S., Gilfanov, M., Churazov, E., & Sunyaev, R. 2006, A&A, 452, 
169 

Revnivtsev, M., Vikhhnin, A., & Sazonov, S. 2007, A&A, 473, 857 
Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2003, A&A, 409, 523 
Schmitt, J. H. M. M., Fleming, T. A., & Giampapa, M. S. 1995, ApJ, 450, 392 
Schmitt, J. H. M. M. 1997, A&A, 318, 215 



- 16 - 

Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Smith, N., & Brooks, K. J. 2008, in Handbook of Star Forming Regions: Vol. II, The 
Southern Sky, ASP Monograph Publications, Vol. 5, ed. B. Reipurth (San Francisco, 
CA: ASP), 138 

Townsley, L. K., et al. 2011, ApJS, submitted (CCCP Intro Paper) 
Townsley, L. K., et al. 2011, ApJS, submitted (CCCP Diffuse Paper) 

Vuong, M. H., Montmerle, T., Grosso, N., Feigelson, E. D., Verstraete, L., & Ozawa, H. 
2003, A&A, 408, 581 

Yonekura, Y., Asayama, S., Kimura, K., Ogawa, H., Kanai, Y., Yamaguchi, N., Barnes, 
P. J., & Fukui, Y. 2005, ApJ, 634, 476 



This preprint was prepared with the A AS macros v5.2. 



-17- 




1.0 



0.8 



_(b)V 



0.6 - 



0.4 



0.2 



0.0 



I'll 

— Simulated foreground stars . 
Simulated background stars 

— Simulated extragalactic 

\ — CCCP sources 
\l 



163 162 161 160 159 163 162 

R.A. (J2000) [deg] 




c 
o 

'■4—' 

D 

4-> 

b 

> 

'4—' 

J2 

E 
u 



160 159 




162.5 



162.0 



161.5 161.0 160.5 

R.A. (J2000) [deg] 



160.0 



159.5 



Fig. 1. — Comparison of histogram (a) and cumulative (b) distributions in right ascension 
among the observed CCCP sources (black), simulated foreground stars (red), background 
stars (green), and extragalactic objects (blue), (c) Spatial distribution of simulated extra- 
galactic objects (blue) superposed on the Chandra Carina exposure map. The yellow con 
tours show the ^^CO emission and green contours indicate C^^O cores from lYonekura et al. 



(120051). 
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Fig. 2. — Comparison of median energy distributions for the observed CCCP sources (black) 
and simulated foreground stars (red), background stars (green), and extragalactic objects 
(blue). 
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Fig. 3. — Panel (a): comparison of J-band distributions for 6182 CCCP sources with avail- 
able 2MASS counterparts (black, magenta, orange, cyan) and simulated foreground stars 
(red), background stars (green), and extragalactic objects (blue). The CCCP sources are 
stratified by median energy: MedE < 1.1 keV (black), 1.1 < MedE < 1.5 keV (magenta), 
1.5 < MedE < 2.5 keV (orange), and MedE > 2.5 keV (cyan). Panel (b): Histograms of 
J-band magnitudes for simulated foreground (red) and ba ckground (green ) stars multiplied 
by the 2MASS J-band detection completeness curve from ISkrutskie et al.l ( l2006l ). As a ref- 
erence, the histogram of the CCCP sources with MedE < 1.1 keV (black) from panel (a) is 
also plotted here. 
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Fig. 4. — Spatial distributions of all observed CCCP sources with 2MASS counterparts (red 
circles) superposed on the Chandra Carina exposure map. Distributions are stratified by 
X-ray median energy, a surrogate for line-of-sight absorption. Panels are stratified by X-ray 
median energies: (a) MedE < 1.1 keV; (b) 1.1 < MedE < 1.5 keV; (c) 1.5 < MedE < 
2.0 keV; (d) 2.0 < MedE < 2.5 keV; (e) 2.5 < MedE < 4.5 keV; and (f) MedE > 4.5 keV. 
The figure legends also provide numbers of plotted CCCP source s. The yellow contours show 
the ^^CO emission and green contours indicate C^^O cores from Yonekura et al. J2005h . 
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Fig. 5. — Same as in Figure HI but for CCCP sources without 2MASS counterparts. 



